Finite-size scaling for the S=l/2 Heisenberg Antiferromagnetic Chain 
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\ Corrections to the asymptotic correlation function in a Heisenberg spin-1/2 antiferromagnetic 

spin chain are known to vanish slowly (logarithmically) as a function of the distance r or the chain 
size L. This leads to significant differences with numerical results. We calculate the sub- leading 
logarithmic corrections to the finite-size correlation function, using renormalization group improved 
perturbation theory, and compare the result with numerical data. 
The correlation function in the spin-1/2 Heisenberg antiferromagnetic chain is difficult to determine from the Bcthc 
[ — ■ Ansatz, so other methods are used for this purpose, such as bosonization and conformal field theory (CFT). These 
methods work very well as a tool to determine long-distance asymptotics. Numerical work however, often 

,— i ■ questioned this approach. This has been explained by the fact that finite size corrections vanish slowly at long 
distances, as 1/lni, and l/ln(r), where L is the system size (periodic boundary conditions assumed) and r the 
S-H ■ separation of the 2 points, due to the presence of a marginally irrelevant operator qQ. 

Following the methods and notation of HH.the continuum limit of the Heisenberg model can be written in SU(2) 
- symmetric form using non-Abelian bosonization Q. The action for the S'C/(2)-symmetric matrix field gj| includes 
the Wess- Zumino term with coefficient k = 1. This theory is equivalent to a free boson defined on a circle. The 
low-energy Hamiltonian in the continuum approximation is given by: 



O 



oo 
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H = H -8it 2 /V3\Jl-Jr, (1) 



O 

where H is the Hamiltonian density for a free boson, 3t,,r are left and right SU(2) currents: 



tr[gTd_ger], J R = —=tr[d + gg 1 er] (2) 



(Here we use the notation A for the marginal coupling constant rather than g as in ||.) 
We now turn to the discussion of the asymptotic correlation function. The spin operators can be written in 
' non-abelian bosonization notation as: 



Sj = (Jl +3r) + const i(-iytr[ S (r} (3) 

Thus the correlation function has uniform and staggered terms, 

G(r) = (S^S*) G u {r) + (-l) r G s {r). (4) 

Both terms vary slowly on the scale of the lattice spacing, and correspond to different Green's functions in the 
Q \ continuum theory In this paper we only consider the staggered term G s (r) , 
O ■ 

^ ; G s (r) (x {tr{* z z){r)tr{<T z z){Q)) L . (5) 

•i-H , 

The staggered correlation function for a finite chain obeys the following renormalization group (RG) equations || : 

[d/dkir + j3(X)d/d\ + 2-y(\)]G s {r,r/L,\) = 0, (6) 
where /3(A) is the beta function for the coupling constant A in Eq.(jl]): 

r/A 



dlnr 



/3(A) (7) 



and 7(A) is the anomalous dimension. In Eq. (^J) the r-derivative acts only on the first argument of G s ; r/L is held 
fixed. Eq. (^) expresses the fact that a rescaling of both lengths L and r by a common factor can be compensated 
for by a change in the effective coupling constant, A(r) and a rescaling of the correlation function. The solution of 
Eq.(|6|) has the form: 
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G s (r, A ) = exp -2 / dlnr' 7 [A(r')] F[r/L, A(r)], (8) 



where A = A(r ) is the "bare" coupling - a coupling at the ultraviolet cutoff scale r , F[r / L, A(r)] is an arbitrary 
function of the effective coupling constant at scale r, A(r). 

The coupling constant flows to zero as the distance r is increased, and one can use perturbative expressions for 
7(A) and /3(A) to determine long-distance properties. The universal terms in the perturbative expansion for the /?- 
function [01 and the anomalous dimension [pUTtl are known, 



/3(A) = -(4tt/V3)A 2 - (l/2)(4^/V3) 2 A 3 (9) 
7 (A) = 1/2 - (tt/\/3)A. (10) 



Thus the effective coupling is given by: 



-L - -L = (An/V3){Hr/r ) + (1/2) ln[ln(r/r )]} + 0(1). (11) 
A(r) A 



Rewriting the integral in Eq.(||) using Eq.(^), one easily finds 

,A(r) 

Thus the Green's function has the expansion 



[ 7 (A)//3(A)]dA = (1/2) In(r/ro) + (1/4) ln[A(r)/A ] + . . . (12) 



G s 



(r,L,A) = i J^ e S: =1 «»[A(rr-AJ] g Fm(r/I)A(r) ™ (13) 



The coefficients, a„ and the functions F m {r/L) can be determined by doing perturbation theory in the bare coupling 
constant and then recasting the resulting expression in terms of the renormalized coupling constant in the form of 
Eq. (|l3|). This automatically incorporates, at low orders, the leading log divergences to all orders in perturbation 
theory. This method is standard in Quantum Chromodynamics calculations and is known as "renormalization group 
improved perturbation theory" . 

The zeroth order term of this sum, Fq(t/L), is given by the free theory - the conformally invariant WZW model 
on a circle of length L, and can be obtained by conformal transformation. Since g has scaling dimension 1/2, for an 
infinite system: 

< tr[^{r)a z {r)]tr[^)(T z ] >= -. (14) 

r 

(Here we have chosen a convenient normalization for the operator g.) For a finite system with periodic boundary 
conditions, we make a conformal transformation from the infinite plane to the cylinder, obtaining: 

< tr[g(r)<T*Mg(0)<r*] > L = , I (15) 

Lsm(7rr/_L) 



Thus we see that: 



Fo(r/L) oc (16) 



Using Eq. (Jll|), we thus obtain the asymptotic correlation function: 



G s{ r, L) - (7r/x)s f n(W£) {Mr/ro) + (1/2) ln[ln ( r/ro)]}. (17) 

Essentially this result was obtained in For related work on the correlation function for the infinite length spin 
chain see Ref. flipf . 

The new result which we derive here is the next order correction, F\(r/L), in Eq. (13). The first order perturbation 
theory result, using the Hamiltonian of Eq. (Ill), gives: 
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S < frfe(r)<7*(r)]trfe(0)<T«] >= | drdxT < tr[g(r, 0)<x z ]ir[g(0, 0)<x 2 ]J L (x, r) • Jb(x,t) > £) (18) 

where T denotes time-ordering. This correlation function can be evaluated using standard CFT techniques. We first 
obtain its value for an infinite system, then obtain the result for finite L by conformal transformation. Using the 
general result for 3-point functions of primary operators we obtain: 

T < tr[g(r, 0K]ir[g(0, 0)<t z ]3 l (x, r) ■ 3 R (x, r) >cx — — — (19) 

(x + zr)(x — r + it)(x — jt)(x — r — «r) 

The normalization constant can be fixed from the operator product expansion (OPE): 

3 l {x,t) ■ J R (x,r)tr[g(0K] -» - J^fffJ,) , (20) 
and the normalization of the zeroth order Green's function in Eq. This gives: 



T < tr[g(r, 0)a z ]tr[g(0, 0)ct z } 3 l (x, t) • 3 R (x, r) > 

r 

167r 2 (x + ir)(x — r + ?'t)(x — ir)(x — r — it) 
The correlation function on the cylinder can be obtained by conformal transformation: 

T < tr[g(r, 0)<r z ]tr[g(0, 0)(r z ]3 L (x, t) ■ 3 R (x, r) > L 

(7r/i) 3 sin(7rr/L) 

167r 2 sin[7r(x + it) / L] sin[7r(x — r + ir) / L] sin[7r(x — it) / L] sin[7r(x — r — it) / L] 
Thus the first order correction is given by: 

5 < ir[g(rK(r)]Mg(0K] >l= (^|) d r £ dx 
(7r/i) 3 sin(7rr/i) 

167T 2 sin[7r(x + it)/L] sin[7r(x — r + it)/ L] sin[7r(x — it)/ L] sin[7r(x — r — it)/ 'L] 



(21) 



(22) 



(23) 



This integral has a logarithmic ultraviolet divergence at r — > 0, and x — > or r. This can be cut off by restricting 
the integral to |x| 2 + r 2 > r\ and |x — r| 2 + t 2 > rg. The resulting logarithmic dependence on tq is exactly what is 
needed in order for the resulting expression to have the form of Eq. (fL3|), since: 



A ° - 1 + ^£ ln(r/ro) . (24 ) 



\jKr) V3 

Remarkably, the integrals can be done exactly. The x- integral can conveniently be done first using contour methods. 
This gives: 



5 < tr[g(7> z (r)]tr[g(0)cr 2 ] > L 

«7T 2 Ao 



2V3L J- 



du 



sinh u sinh(u + iirr/L) 



complex conjugate 



(25) 



(Here u = 2ttt/L and the integral is cut off at \u\ > r^/L.) This indefinite integral may be done exactly by changing 
variables to tanhu, giving: 

<tr[g{r)<T z {r)}tr{g(0)cT z } >L 

' -{1 + ^^{ln[(i/r ) sin(7rr/L)] + constant}}. (26) 



(L/tt) sin[7rr/i] ^3 

This is the result of perturbation theory to first order in the bare coupling constant Ao. The next step is to "renor- 
malization group improve" this result by matching it to the expression in Eq. (IK})) . Expanding this expression to first 
order in the bare coupling constant, using Eq. (pdj) gives: 
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G s (r) oc 



1 



(L/tt) sin[7rr/L] 
Comparing Eq. (p6h and (Et]) we see that: 



[1 + A (27t/V3) In(r/r )][l + \ F 1 (r/L)/F (r/L)}. 



2n 



F 1 (r/L)/F Q (r/L) = ^= In 



~L . 


/"rrr\ 


— sm 




r 


<~LJ 



constant. 



Thus our renormalization group improved expression for the correlation function is: 



G s (r,L,X ) oc 



1 



An 



{L/tt) sin[7rr/L] y A(r) 



L ( Ttr\ 

- sm {T) 



constant 



(27) 



(28) 



(29) 



The advantage of this RG improved expression is that we may now go to arbitarily large r, a limit in which the large 
logarithms, ln(r/ro) invalidate finite-order perturbation theory and infinite resummations of most divergent diagrams 
are neccessary. This is automatically taken care of by Eq. (|29j) together with the expression for A(r) in Eq. (y). In 
this asymptotic limit we may use: 



A(r) 



V3 



4-7T ln(r/ro) ' 

for the factor of A(r) inside the curly brackets in Eq. (|29[). This gives: 
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(30) 



Gs(r) 



(L/tt) sin[7rr/L] 
1 1 



[ln(Cr/r ) + (l/2)ln[ln(r/r )]] 1 / 2 



21n(r/r ) 



In 



L ( irr 
— sm — 
r V L 



+ constant 



where: 



C 



a \/3/47rAo+0(l) 



(31) 



(32) 



We now see that the "constant" term can be adsorbed, to lowest order in l/ln(r/Vo) , into a rescaling of C (i.e. a 
shift of An) so we henceforth drop it. This is all information that can be extracted from the RG to this order. Two 
non-universal free parameters remain: the overall amplitude, A and the bare coupling An, appearing as the constant 
C in Eq. (j3^). The amplitude A was recently determined, from Bethe ansatz results [[[o) for the S=l/2 Heisenberg 
model, to be @ A = (2tt)~ 3 / 2 . The bare coupling Ao (for some conveniently chosen value of ro) is not known exactly 
for the S=l/2 Heisenberg model. It can be determined by fitting numerical results. Since the same bare coupling 
constant also appears in various other finite size corrections, including the finite size spectrum |(| various consistency 
checks could be made. However, this requires further calculations to ensure consistency of the cut-off schemes in 
various calculations and we do not attempt it here. 

Note that if we chose z = (L/tt) sm(irr/L) instead of r as the length scale in Eq.(||), we would find, to this order, 
that the finite-size correlation function at distance r is given by an infinite chain correlation function at z: 



G s (r,r/L) = G s (z,0) 



r(27r) 3 / 2 



{ln(Cz/r ) + (l/2)ln[ln(2/r )]} 1/2 . 



(33) 



However, this fact does not necessarily hold for higher order RG. As follows from our result Eq.(|3l|), the finite-size 
correction to the asymptotic correlation function vanishes as r 2 / L 2 as one approaches the infinite-chain limit, r/L — > 0. 
The expansion in l/ln(r/rn) remains, but the functions F m (r/L) (see Eq.(|l3|) approach constants in this limit. (We 
have only checked this for Fq(t/L) and Fi(r/L)) 

It is interesting to compare our result Eq.(31) with phenomcnological expressions used to fit numerical data. Koma 
and Mizukoshi |5| used the scaling function of the form: 



G,(r,L) = 



A{ln[(L/7rr )sin( 7 rr/L)]} 1 / 2 
{L/tt) sin(7rr/ ' L) 



(34) 



with A — (27r)~ 3 / 2 ~ 0.0635, close to the exact answer. This form is equivalent to Eq.(|33]) with the slow log log 
term replaced by a constant. The best fit was obtained for A ~ 0.065, close to the exact answer. Kubo et al |l| and 
Hallberg et al [0 defined a scaling function f(r/L), 
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G a (r,L) = G a (r,oo)f(r/L), 
and adopted a phenomenological expression, 

f(x) cx [l + J 4sinh 2 (Ba;)] 2 ' ) . 



(35) 



(36) 



Note that Eq. (poj) does not agree with our RG analysis. They found the best fit for A = 0.28822, B = 1.673, 
2rj = 1.805. As it has been noted by one of us, the form is remarkably (within 0.057.) close to the CFT prediction for 
the general xxz model, since 



1 + 0.28822 sinh 2 (1.673a;) 



1/2 



sin(7nr) 

Further, for 77 close to 1 we can replace the phenomenological formula Eq.(^6|) by an equivalent, 



[1 + 0.28822 sinh 2 (1.673a;)] 



2»J 



sin(7ra;) 



1 + (77-I) In 



sin(7rx) 



The scaling function f{x) is similar to F(X(r),x) defined in Eq.(|3l"|). We find that rj depends on r, 

1 



7](r) 



1 



21n(cr) 



(37) 



(38) 



(39) 



Finally, let us compare our theoretical expression Eq.(31) with DMRG data of Ref. M. We have found that the 
log log term, which is higher order, is almost constant, and does not influence the comparison. Since adding a loglog 
term introduces one more free parameter, we decided to drop it. To obtain data collapse for G s (r,L) in a one-loop 
RG we use Eq. (|33|) , with the log log term dropped and use C as a free parameter. We compare this with a zero-order 
non-interating result, 



Gf ree (z) — const 



1 



(40) 



The result is shown in Fig.JT]. Both expressions use one free parameter. The better fit produced by our one-loop 
expression is obvious. 




FIG. 1. Scaled DMRG data for spin correlation function compared with one-loop theoretical expression (solid line) and free 
boson result (dashed line) 

Alternatively, we can use the result Eq.(|3l|), with A = l/(27r) 3 / 2 , the log log term dropped, and C taken as a free 
parameter. (These two expressions are the same to the order that we have calculated in 1/lnr.) The comparison of 
this formula with numerical data is shown in Fig.^. For comparison we also show the result without the correction 
that we have calculated, 



G s (r) 



AMCr/ro)] 1 ' 2 
(L/ir) s'm[irr/L] 



(41) 
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FIG. 2. Numerical data for the spin-spin correlation function G s (r,L) from Ref. |4J versus r/L for r = 7 — 27. Solid lines - 
our one- loop result (with one parameter - c), dashed line - the result without our correction (also with one parameter - c). 
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FIG. 3. Spin-spin correlator multiplied by r, rG s (r, oo), vs r extracted from the DMRG data of Ref. Solid line - the 
one- loop result. 



It is instructive to find the value of A from the numerical data to compare it with exact value. Using Eq.(|3I 
with A and C as free parameters to fit the numerical data, we find A = 0.0636427, which differs only 0.2% from the 
exact answer. This is much better than what we would have found without the correction calculated in this paper, 
A = 0.0578896, which is off by some 97„. 

Finally, we can extract from the data the behavior of the correlation function at L — oo, using Eq.(^) and dividing 
it by 



F (r/L) + Fi(r/L)A(r) = 1 + 



2 ln(cr 



■In 



L /7tT\ 

— sin — 

7rr V L / 



L sin[7rr/L] 



The result is given by 



G s (r, oo) = 



(2tt) 3 /V 



(42) 



(43) 
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with the same constant c. This way of extracting L — oo behavior from finite-size data is different from the phe- 
nomenological scaling expression that Hallberg et al. used in Ref. Q, G s (r, L) — G s (r, oo)f(r/L). The result is shown 
mFigf 

In conclusion, we have shown that bosonization approach provides an accurate description of the spin correlation 
function in a finite spin-1/2 Heisenberg spin chain. Our theoretical result Eq. ([?!]) compares favorably with numerical 
data at long length scales. We are greatful to K. Hallberg for providing us numerical data from Ref. This work 
was supported by NSERC of Canada. One of us (VB) acknowledges support by the NHMFL through NSF cooperative 
agreement No. DMR-9527035 and the State of Florida. 
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